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POPULAR SUMMARY 

Principle Component Analysis of the Evolution of the Saharan Air Layer and Dust 
Transport: Comparisons between a Model Simulation and MODIS Retrievals 


Every summer, large quantities of mineral dust are blown westward over the tropical 
Atlantic Ocean from the Saharan desert. The frequency of these dust outbreaks is 
modulated by the occurrence of African Easterly Waves, which propagate from the west 
coast of Africa across the tropical Atlantic Ocean. The dust is transported in an elevated 
(2-5 km altitude) layer of warm, dry air called the Saharan Air Layer (SAL) that moves 
above the relatively cool, moist marine boundary layer over the ocean surface. Recent 
investigations have li nk ed the transport of the SAL with suppressed tropical cyclone 
activity in this region: the warm air anomalies within the SAL serve to stabilize the 
atmosphere, suppressing the convective activity which drives the tropical cyclone activity 
and plays a role in, for example, the development of hurricanes. Other investigations 
have demonstrated a spatial correlation of the warm air anomalies in the SAL with 
enhanced dust concentrations. In short, the radiative heating within the elevated dust 
layer plays a role in suppressing tropical cyclone activity. Numerical models (e.g., 
forecast models) that do not appropriately simulate the development of the SAL will have 
a tendency to overestimate the intensity of associated tropical cyclones. 

In this study we investigate the evolution of the SAL and associated dust transport during 
the summer of 2002 using a combination of meteorological datasets, numerical models, 
and satellite observations. We apply the statistical technique of principle component 
analysis to meteorological data from the National Center for Evironmental 
Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalyses in order 
to diagnose the onset and direction of transport of SAL warm anomalies. Using the 
NASA/University of Colorado Community Aerosol and Radiation Model for 
Atmospheres (CARMA) and the NCEP/NCAR reanalyses we simulated the dust 
transport over this region for the same period. The validity of the model was assured 
through comparisons with dust aerosol observations made with the NASA Moderate 
Resolution Imaging Spectroradiometer (MODIS) onboard the polar-orbiting Terra 
satellite. We find that near the west African coast the presence of dust is spatially well 
correlated with the warm anomalies in the meteorological datasets, but that this 
correlation breaks down further westward. We also diagnose two main transport 
pathways for the dust, one directly west of Africa, and a second more northwestward, 
toward the mid-latitude Atlantic and Europe. The statistical analysis of the warm 
anomalies provides an index for diagnosing which pathway is more important at a given 
time. We further find a significant linkage between the evolution of the SAL and the 
Western African Monsoon, as diagnosed by the difference between high altitude winds 
and the winds in the SAL, although the physical mechanism for this correlation remains a 
subject of ongoing investigation. Furthermore, the apparent breakdown in the spatial 
correlation between warm anomalies and dust concentrations in the western tropical 
Atlantic needs to be investigated more. Although this discrepancy may be subject to 
dynamical controls such as the northward advection of warm air parcels from the deeper 



tropical atmosphere, it may also be somewhat mitigated by the radiative heating of the 
dust, which is not accounted for in the current model and may not be properly accounted 
for in the meteorological analyses because of sparse data in this region, but can be 
considered in the context of newer NASA models for dust transport and data assimilation. 
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Abstract 


The onset and evolution of Saharan Air Layer (SAL) episodes during June-September 
2002 are diagnosed by applying principal component analysis to the NCEP reanalysis 
temperature anomalies at 850 hPa, where the largest SAL-induced temperature anomalies 
are located. The first principal component (PC) represents the onset of SAL episodes, 
which are associated with large warm anomalies located at the west coast of Africa. The 
second PC represents two opposite phases of the evolution of the SAL. The positive 
phase of the second PC corresponds to the southwestward extension of the warm 
anomalies into the tropical-subtropical North Atlantic Ocean, and the negative phase 
corresponds to the northwestward extension into the subtropical to mid-latitude North 
Atlantic Ocean and the southwest Europe. A dust transport model (CARMA) and the 
MODIS retrievals are used to study the associated effects on dust distribution and 
deposition. The positive (negative) phase of the second PC corresponds to a 
strengthening (weakening) of the offshore flows in the lower troposphere around 10°- 
20°N, causing more (less) dust being transported along the tropical to subtropical North 
Atlantic Ocean. The variation of the offshore flow indicates that the subseasonal 
variation of African Easterly Jet is associated with the evolution of the SAL. Significant 
correlation is found between the second PC time series and the daily West African 
monsoon index, implying a dynamical linkage between West African monsoon and the 
evolution of the SAL and Saharan dust transport. 
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Introduction 

The Saharan Air Layer (SAL) is a well-mixed layer of dry, warm air located above the 
marine boundary layer with the top extending to ~500 hPa [Carlson and Prospero, 1972], 
From late spring to early fall, it is advected westward with Saharan dust across the North 
Atlantic Ocean [Carlson and Prospero, 1972; Prospero and Carlson, 1972]. Recently, 
Dunion and Velden [2004] used the Geostationary Operational Environmental Satellite to 
track the SAL across the North Atlantic Ocean. They found a phenomenological 
connection between the SAL and the suppression of the tropical cyclone (TC) activity. 
They indicated that the dissipation of deep convection occurring in the interior of the 
SAL plays a role in suppressing the development of TC activity. Using satellite-derived 
aerosol optical thickness and cloud brightness temperature, Wong and Dessler [2005] 
verified the suppression of deep convection associated with the warm and dry anomalies 
of the SAL. A model that does not appropriately simulate the SAL will therefore tend to 
overestimate the intensity of TCs that are interacting with the SAL [Dunion and Velden, 
2004]. Understanding the properties of the SAL and studying its evolution are therefore 
important tasks for numerical prediction of TCs. 

The SAL contains Saharan dust, which propagates with the SAL westward across the 
North Atlantic Ocean. There have been many model studies about large-scale 
intercontinental transport of dust. Tegen and Miller [1998], Mahowald et al. [2003], and 
Colarco et al. [2003] all demonstrated that the downwind dust distribution is mainly 
influenced by the dynamical transport rather than the source strength. 

The relationship between the dynamical and dust fields has been studied previously. 
Alpert et al. [1998] showed that there is a spatial correlation between the reanalysis warm 
increments and the dust concentrations in the eastern tropical North Atlantic Ocean. A 
model-simulated climatology by Jones et al. [2004] indicated that the westward transport 
of dust is modulated by African Easterly Waves, which propagate from the west coast of 
Africa across the tropical Atlantic Ocean. Karyampudi et al. [1999] used lidar 
observations to show detailed spatial relationships between the dust plumes and the 
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dynamical fields, with the maximum dust concentration located close to the mid-level jet. 
These studies all imply a coherency between the migration of the SAL and the Saharan 
dust, and that the linkage between them is related to the propagation of African Easterly 
Waves. 

With the availability of data from the Moderate Resolution Imaging Spectroradiometer 
(MODIS) [Kaufman et al., 1998; King et al., 2003], sophisticated statistical analyses can 
be performed on global-scale satellite data to reveal new information about the 
interaction between the large-scale dynamical fields and intercontinental dust transport. 
Kaufman et al. [2005] applied the aerosol optical thickness (AOT) from MODIS to 
estimate the cross-Atlantic flux and deposition budgets of the Saharan dust. They 
showed a high correlation during the summertime between the monthly averaged MODIS 
dust AOT and westward wind speed at 750 hPa around the tropical-subtropical west coast 
of Africa, verifying that the dust outflow from the African continent to the eastern North 
Atlantic Ocean is modulated by the offshore easterlies during summer. 

In this work, we study the evolution of the SAlL and the associated dust transport and 
deposition over subseasonal time scales using both model simulations and measurements 
from the MODIS. In Section 2, we describe the dust transport model and the 
meteorological data used to drive the model. We also describe the satellite retrievals 
being used to evaluate the model results. In Section 3, we apply statistical techniques to 
define two indices describing the onset and the evolution of each SAL episode during 
June-September 2002. We identify different features of the evolution of the SAL and the 
related pattern of West African circulation. We also estimate the relative distribution of 
dust deposition between the eastern tropical North Atlantic and the northwest Africa to 
southwest Europe. Conclusions are presented in the last section. 

Model and Data 


1. Model Simulation 



5 


Global dust distributions are simulated for the period June-September 2002 using the 
offline Community Aerosol and Radiation Model for Atmospheres (CARMA, Toon et al. 
[1988], Ackerman et al. [1995]), following a methodology similar to that in Colarco et al. 
[2003]. Briefly, the NCEP/NCAR reanalysis package [Kalnay et al., 1996] and the 
NCAR Model for Atmospheric Transport and Chemistry (MATCH, Rasch et al. [1997]) 
provide the meteorology and physical parameterizations (e.g., clouds, convection) 
necessary to simulate the atmospheric state. Dust emissions are based on soil credibility, 
soil wetness, and 10-m wind speeds [Ginoux et al. 2001]. Dry deposition is the most 
important dust removal mechanism in terms of the total aerosol mass, and is highly 
dependent on particle size, so to accurately simulate the evolution of the dust particle size 
distribution we distribute the dust over eight size bins between 0.1 and 10 pm radius and 
compute the fractional emissions for each bin. Dust emissions and meteorology are 
archived every six hours over the simulation period. The horizontal spatial resolution of 
these fields and the subsequent dust simulations is approximately 1.875° x 1.875° on a 
regular longitude and Gaussian latitude grid. The vertical grid is resolved in 28 sigma 
layers extending from the surface to approximately 35 km altitude; 18 of these layers are 
below 10 km altitude. These fields are read into CARMA, which merges the emissions 
and meteorology with the aerosol microphysical and transport continuity equation to 
simulate the full dust lifecycle. The timestep in CARMA is 1800 seconds, so input fields 
are linearly interpolated to the model current time. Dust removal is by dry deposition, in- 
and below-cloud wet scavenging in large-scale cloud systems, and wet scavenging in 
convective updrafts [Balkanski et al., 1993]. The dust AOT at 550 nm is computed for 
each model column from the CARMA column-integrated particle size distribution 
assuming Mie scattering and a refractive index of 1.50 - 0.01/. Further details of the 
model and its validation are presented in Colarco et al. [2003]. 

2. Satellite Data 

The MODIS aerosol algorithm inverts the measured radiance from six MODIS bands 
(0.55-2.1 pm) to retrieve aerosol properties at 10-km resolution [Kaufman et al., 1997; 
Levy et al., 2003; Remer et al., 2002]. Radiance at these wavelengths comes almost 
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entirely from reflected solar radiation; as a result, these measurements are only available 
during daytime. Daily statistics of aerosol properties at grids of l°xl° resolution are 
archived in the MOD08 level-3 product (version 4) [King et al., 2003]. We analyze here 
the 0.55-pm AOT product measured between June and September 2002. We use 
measurements made by the MODIS onboard the Terra satellite, which crosses the equator 
around just before local noon. In this paper, we focus on the comparison over oceans 
only, where the retrieved AOTs are more reliable than those over lands. 

MODIS AOT include radiance extinction from all kinds of aerosols. In the tropical and 
subtropical North Atlantic Ocean, the main contributions of AOT are from dust, maritime 
aerosol, and aerosol from biomass burning in South Africa. Kaufman et al. [2005] 
estimated the seasonal cycle of the relative contributions from these three types of 
aerosols to the total AOT. In the region of tropical to subtropical North Atlantic Ocean, 
contribution from dust to the total AOT is more than 50% during summer to early autumn 
(June-September). The remaining contributions are mainly from maritime aerosol. Since 
the CARMA simulation discussed in this paper includes only dust, we would expect a 
difference of about a factor of 2 between the MODIS and CARMA AOT. The actual 
difference is larger because CARMA in general underestimates AOT. In this study, 
however, the important feature is the patterns of dust distribution. Differences between 
MODIS and CARMA AOT values do not affect our main conclusions. 

The warm anomalies associated with the SAL obtained from NCEP temperatures can be 
compared by those obtained from the Atmospheric Infrared Sounder (AIRS) [Strow et al., 
2003; Susskind et al., 1998] onboard the Aqua satellite. AIRS measures upwelling 
infrared radiance in 2378 spectral channels and uses those measurements to provide 
accurate vertical profiles of atmospheric temperatures at a horizontal resolution of ~50- 
km at nadir. The temperature retrievals are validated to 1 K in 1-km vertical layers over 
ocean [http://www-airs.jpl.nasa.gov/press/news_factsheets.html]. In this study, we used 
the daily level-3 product (version 3) with l°xl° spatial resolution to evaluate a SAL 
episode that occurred in early September 2002. 
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Results 

1. Principal Components and Time Evolution 

We analyze the NCEP temperature data at 850 hPa, where the warm anomalies associated 
with the SAL reach their maxima, in order to investigate the evolution of the SAL. We 
first remove the seasonal cycle of the temperatures at 850 hPa for the summertime of 
2002 (June 1 to September 30). The seasonal value of temperature at each grid on a 
given day is estimated by taking the 31 -day average of the temperatures from the -15 to 
15 days counting from the given day. The anomalies are the difference between the 
noon-time temperatures and the seasonal value for each day and each grid. Positive 
(negative) anomalies mean that the temperatures are higher (lower) than the seasonal 
values. We then apply principal component analysis (PCA, [Bretherton et al., 1992; 
North et al., 1982; Wong and Wang, 2003]) to the temperature anomalies in the region of 
0°-80°W and 0°-40°N, which covers the areas of tropical to mid-latitude North Atlantic 
Ocean and the west coast of Africa. The time series and the patterns of the empirical 
orthogonal functions (EOFs) for the first two principal components (PCs) are shown in 
Fig. 1 and the left two panels of Fig. 2. The first PC represents about 18% of the 
variability of temperature anomalies over the region, and the second PC represents about 
11%. The errors of the eigenvalues are estimated by the rule stated by North et al. [1982], 
The eigenvalues and the associated errors for the first, second, and third PCs are 
6.67±0.85, 3.82±0.49, and 3.12±0.40, respectively. Although the spread of the second 
and third PCs are marginally overlapped, the comparisons shown below between satellite 
retrievals and CARMA dust simulation justify the physical reality of the second PC. 
Similar analysis using summertime NCEP temperatures for longer period (1991-2002, 
not shown) also gives similar pattern for the second PC. 

We also analyze the NCEP wind fields and the distributions of dust AOT from the model 
associated with the two PCs. The altitude at which the winds are most correlated with the 
transport of Saharan dust out of the African coast is about 700 hPa [Kaufman et al., 2005]; 
therefore, we analyze the NCEP winds at 700 hPa. We estimate the seasonal values of 
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the winds (and dust AOT) at each grid on a given day by the 31 -day averages of 
northward and eastward components (and dust AOT) from -15 to 15 days counting from 
the given day, similar to what we do for the temperatures. The wind anomalies (hereafter, 
referred to as the transient eddy winds) and the dust AOT anomalies are computed as the 
residual components after the removal of the seasonal values. The wind anomalies reflect 
the subseasonal variations from the seasonal values. Positive (negative) AOT anomalies 
mean that the concentrations of dust are higher (lower) than the seasonal values. The 
dust AOT anomalies and the eddy winds are then regressed on the time series of the two 
PCs, in order to illustrate how the dust concentrations and the winds are associated with 
the temperature anomalies. The results are shown in the right two panels of Fig. 2. 

The first EOF (EOF-1, Fig. 2 top left) shows large warm anomalies located off the west 
coast of Africa. The warm anomalies are centered at about 30°N and extend toward 
Europe as well as into the tropical-subtropical North Atlantic Ocean. The regression of 
the dust AOT anomalies on the first PC (Fig. 2 top right) is spatially well correlated with 
the pattern of EOF-1. However, the maximum of dust AOT anomalies is located slightly 
to the south of the warm anomalies and coincides with the location of the strong off-shore 
eddy winds to the east and south of the maximum warm anomalies. This is in agreement 
with the lidar measurements showing that dust anomalies tend to be close to the location 
of mid-level easterly jet [Karyampudi et al., 1999]. In the following, we will utilize the 
time series of the first PC (Fig. 1 top) as an index to indicate the timing of the beginning 
of SAL episodes at the west coast of Africa. We define the onset of SAL episodes as 
when the index exceeds the value of one (meaning that the warm anomalies are more 
than the lo-confidence level). 

The second EOF (EOF-2, Fig. 2 bottom left) shows a dipole pattern. One end of the 
dipole (positive values) is located at the west coast of Africa between about 15°N and 
30°N, at the south of the maximum seen in EOF-1. The other end of the dipole (negative 
values) is located northward of 30°N over the northwest- African coast and southwest 
Europe. The regression of dust AOT anomalies on the second PC also shows a similar 
dipole pattern (Fig. 2 bottom right) with the centers located slightly southward of the 
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centers of the temperature anomalies. The second PC represents two different directions 
for SAL migrations during June-September 2002. We utilize the time series of the 
second PC (Fig. 1 bottom) as an index to define the opposite phases of the evolution of 
the SAL. The positive phase is when the index has positive values, corresponding to the 
southward and westward migration of the SAL into the tropical-subtropical North 
Atlantic Ocean; the negative phase is when the index has negative values, corresponding 
to the westward and northward migration of the SAL into the subtropical to mid-latitude 
North Atlantic Ocean and the southwest Europe. 

The SAL’s two different migration pathways are related to the eddy wind pattern off the 
west coast of Africa. For the positive phase (Fig. 2 bottom right), the northeasterly eddy 
winds at the east and south boundaries of the warm anomalies around 15-20°N strengthen 
the easterly trade winds transporting dust westward. For the negative phase, the direction 
of eddy winds is reversed and becomes blowing onshore, weakening the offshore 
transport of dust at 15-20°N into the Atlantic Ocean. Instead, more dust tends to be 
distributed northward to northwest Africa and southwest Europe. In conclusion, the 
migration of SAL from the African continent is correlated with the variation of the 
offshore winds at about 15-20°N over the west coast of Africa. Previous studies also 
found high correlations between dust AOT and westward winds off the west coast of 
Africa at the same latitudes [Karyampudi and Carlson, 1988; Kaufman et al., 2005; ]. 

2. Temporal Evolution 

We pick from Fig. 1 four individual SAL outbreaks, starting on June 10 (I), June 21 (II), 
July 23 (III), and September 5 (IV). All four episodes have the values of the first PC 
index greater than one (meaning that the warm anomalies are more than l-o confidence 
level at the onset). We examine the evolution of each episode by computing the time-lag 
composites of the temperature and dust AOT anomalies. The time-lag composites are 4- 
day averages of the anomalies with time lags from the onset dates. For example, for the 
episode (IV) the zero day-lag composite is the 4-day averages of the anomalies from 
September 5 to 8, the one day-lag composite is the 4-day averages from September 6 to 9, 
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the two day-lag composite from September 7 to 10, and so on. The composites of the 
model variables allow us to compare the model results with the satellite data, because 
missing values in the satellite observations due to the clouds or the orbital gaps are 
mostly mended by the averages. 


Since the temperature anomalies on each day can be expanded by the PCs, we can use the 
time series of the first two PCs to give an overview of the evolution of the warm 
anomalies near the west coast of Africa. Let T(x,t i ) be the temperature anomalies at 

position x on day t , , the 4-day averages of the anomalies starting on the day t i can be 
expressed as 


T (x, ^ ) = 7 Z T (x, ',+«) = 7 Z Z c « + »)•</„(*) = Z f 7 Z c m (h +n )\v m 0) 


n = 0 


m «=0 


where c m (t) represents the time series for the m-th PC and i// m (x) represents the m-th 


l 3 

EOF. Therefore, the time-lag composites of the PC time series, — Z c mft + n ) > provide a 

4 n = 0 

means to investigate the evolution of each episode. Figure 3 shows the time-lag 

l 3 

composites of the first two PC time series (i.e., — Z c m ft + n )> for m = 1,2 ), beginning 

4 n = 0 


on the individual onset dates of the four main episodes. Episode I has relatively flat time- 
lag variation of the composites of the first PC time series, because of the several 
successive events after the first episode (e.g., events from June 13 to 15, see Fig. 1 top 
panel). The time-lag composites of the second PC time series are mainly negative for 
Episode I, meaning that the eddy winds near the west coast of Africa are blowing onshore, 
and are mainly positive for Episode II, meaning that the eddy winds are blowing offshore. 


Figure 4 shows the corresponding spatial distributions of the temperature anomalies and 
CARMA dust AOT anomalies for the second-day lag of Episodes I and n. The 
distributions of the warm anomalies of these two episodes (Fig. 4 left) are consistent with 
the EOF-2 pattern (Fig 2. bottom left). The negative-phase dominated Episode I has 
warm anomalies mainly distributed north of 20°N over the subtropical to mid-latitude 
Atlantic Ocean. The positive-phase dominated Episode II has warm anomalies extending 
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southwestward into the tropical Atlantic Ocean. In Episode I, the transport of dust along 
15-20°N is suppressed by the onshore eddy winds at the west coast of Africa. 
Consequently, more dust tends to travel northward to Europe, and some heads 
southwestward into the deep tropical Atlantic Ocean (Fig. 4, top right). The anomalies 
seen in the central Atlantic Ocean around 15°-20°N are the reminiscent from a previous 
episode. In contrast, during Episode II, the dust is transported into the Atlantic Ocean 
around 15-30°N by the offshore eddy winds (Fig. 4, bottom right). 

For Episode III, the time-lag composite of the second PC time series (Fig. 3 bottom left) 
changes sign around the second-day lag, implying that there is a change in the near-coast 
wind direction, accordingly altering the evolution of the warm anomaly. For Episode IV, 
the 4-day averaged second-PC time series remains positive for 4 to 5 days of time lag. 
For this period, the offshore eddy winds at the west coast of Africa dominate for the 
evolution of the SAL. On the 5-th to 6-th day of time lag, the phase becomes negative 
and onshore eddy winds dominate. In the next three subsections, we will focus on 
Episodes III and IV and compare the reanalysis and model results with satellite retrievals. 

3. Evolution of Episode III 

Figure 5 illustrates the evolution of the warm anomalies associated with the SAL for 
Episode III (onset day on July 23). The spatial distributions of the warm anomalies are 
qualitatively similar to the pattern of the SAL tracking imagery shown by Dunion and 
Velden [2004] (comparing the Lag 2 and Lag 4 figures with their Fig. 1), although they 
showed the imagery in a different year. Figure 5 also shows the eddy wind field 
composites associated with the warm anomalies on each day lag. The corresponding 
CARMA and MODIS dust AOT composites are shown in Figs. 6 and 7. 

The composite on the 0-day lag (July 23-26) shows that the SAL warm anomalies at 850 
hPa extend southwest toward the tropical Atlantic Ocean (Fig. 5 top left). The 
corresponding composite of the second PC time series (Fig. 3 bottom left) is positive and 
the distribution of the warm anomalies is consistent with the positive-phase of the second 



12 


PC. Over the eastern Atlantic Ocean near the west coast of Africa, offshore eddy winds 
located at the south of the warm anomalies are responsible for blowing the dust off the 
west coast of Africa; hence, the dust is located slightly south of the warm anomalies. 
However, the simulated dust AOT anomalies (Fig. 6 top left) extend more westward than 
the warm anomalies. The MODIS AOT measurements (Fig. 7 top left) confirm the 
distribution of the CARMA dust AOT, although the simulated dust AOT anomalies are 
lower than the MODIS retrieved values. The incoherence in the spatial distribution of the 
warm anomalies and the dust AOT anomalies in the tropical North Atlantic may be 
attributed to the pattern of the eddy winds. Dust AOT anomalies tend to follow the 
westward winds, while the warm anomalies tend to follow the northward winds that bring 
in warm air parcels from the south. Therefore, these will generally not be collocated, 
implying that dust is not a good proxy for the SAL in the regions far away from the 
African continent [Dunion and Velden, 2004; Wong and Dessler, 2005]. 

The evolution of the SAL and the associated dust transport undergo changes drastically 
as the phase of the second PC shifts from positive to negative (after the second-day lag in 
Fig. 3, bottom left). The warm anomalies on the west coast of Africa are pushed 
northward by a strengthening cyclonic circulation located south of the warm anomalies 
near the west coast of Africa around 20°N (Fig. 5 top right). The transport of dust is 
accordingly altered by the changed circulation. On the fourth-day lag, the eddy winds at 
the west coast of Africa around 20°N blow onshore, and the dust follows the circulation 
pattern and is distributed northward over the northwest coast of Africa into the mid- 
latitude Atlantic Ocean (Figs. 5 and 6 bottom left). The corresponding MODIS AOT 
images (Fig. 7 top right and bottom left) also capture the northward transport of dust. On 
the sixth-day lag, the warm anomalies (Fig. 5 bottom right) and the dust (Fig. 6 bottom 
right) are located near the northwest African coast as is the MODIS AOT image (Fig. 7 
bottom right). 


4. Evolution of Episode IV 
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Figure 8 shows the evolution of the warm anomalies of Episode-IV (onset on September 
5) and the associated eddy wind composites. We also compare in Fig. 9 the evolution of 
this episode with that seen in AIRS retrieved temperatures for the second and fourth-day 
lags. Since we have only AIRS data for September of 2002, we obtain the warm 
anomalies in Fig. 9 by subtracting the monthly mean temperatures of September 2002, 
instead of the averages of the -15 to +15 days from the onset day. Comparing the lag-2 
and lag-4 panels of Fig. 9 with the corresponding panels of Fig. 8, we can see that the 
various methods of determining the anomaly field produce quantitatively similar results. 
Fig. 9 shows that the evolution of the warm anomalies captured by the NCEP reanalysis 
temperatures is remarkably consistent with that by the AIRS retrieved temperatures. In 
this episode, the warm anomalies move westward along the latitude of about 20°-30°N. 

The associated CARMA and MODIS dust AOT anomalies are shown in Figs. 10 and 11, 
respectively. On the 0- and second-day lags, the SAL warm anomalies are located just 
off the west coast of Africa. The offshore eddy winds at the south boundary of the warm 
anomalies (Fig. 7 top) are again responsible for blowing the dust off the west coast of 
Africa, and the maxima of the dust AOT anomalies are located slightly southward of the 
maxima of warm anomalies. On the second- to fourth-day lags, the phase of the second 
PC is still positive (Fig. 3 bottom right), and the warm anomalies and the dust are moving 
westward into the tropical-subtropical North Atlantic Ocean by the enhanced easterly 
trade winds. 

On the fourth- and sixth-day lags, the warm anomalies have detached from the west coast 
of Africa and been moving toward the central Atlantic Ocean (Fig. 8 bottom panels). 
This feature is captured by the AIRS measurements (Fig. 9 bottom panels). The 
corresponding westward transport of dust is captured by the CARMA simulation (Fig. 10 
bottom panels) and the MODIS measurements (Fig. 1 1 bottom panels), although detailed 
difference in dust distributions between CARMA and MODIS becomes more evident in 
later time-lag composites. 


5. Distribution of Dust Deposition 
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According to the above discussion, the distribution of the deposition of dust over the 
eastern tropical North Atlantic Ocean and the northwest Africa to southwest Europe is 
highly dependent on the phase of the second PC. To quantify this, we compute the 
distribution of dust deposition between two regions using the CARMA simulated dust 
deposition, which is the sum of dry deposition and scavenging in large-scale and 
convective clouds. The first region covers the eastern tropical North Atlantic Ocean 
(20°-40°W and 0°-30°N). The second region covers the northwest Africa to southwest 
Europe (0°-20°W and 30°-45°N). 

Table 1 shows the relative distribution of dust deposition between the two regions. 
Because of the change of phase in the second PC during Episode III, we average the dust 
deposition for the first four days (day-0 to 3) and the successive four days (day-4 to 7) 
from the onset dates of each episode. The relative distributions are shown in percentage 
of the combined dust deposition in the two regions for each averaged period in each 
episode. Also shown in Table 1 are the combined dust depositions of the two regions and 
the averaged second PC values for the averaging period. 

In Episode I, the second PC is mainly in the negative phase and relatively more dust is 
deposited into northwest Africa and southwest Europe. Episode II is mainly in the 
positive phase of the second PC and has more dust deposited over the eastern tropical 
Atlantic Ocean. In Episode III, the phase of the second PC changes from positive to 
negative around the second-day lag (Fig. 3 bottom left). The first four-day average for 
Episode III in Table 1 shows dust deposition distributed mainly over the eastern tropical 
Atlantic Ocean (85% vs. 15%); however, the second four-day average shows an increased 
portion of dust deposition distributed over the northwest Africa and southwest Europe 
(50% vs. 50%). The second PC for Episode IV does not change its sign until the fourth- 
day lag (Fig. 3 bottom right), when most dust associated with the warm anomalies has 
already migrated into the central North Atlantic Ocean. Therefore, a large portion of dust 
deposition is distributed over the eastern tropical Atlantic Ocean. 
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In conclusion, there are two transport pathways for North African dust. The first is 
directly westward into the tropical to subtropical Atlantic, while the second is 
northwestward toward mid-latitude Atlantic and Europe. The second PC of the 
temperature anomaly gives us a measure of the importance at any given time of this 
pathway: more (less) dust is deposited over northwest Africa and southwest Europe than 
over the eastern tropical North Atlantic Ocean, if the second PC has negative (positive) 
phase. 

Conclusions and Discussion 

Principal component (PC) analysis is applied to NCEP temperature anomalies at 850 hPa 
for June-September 2002. The first PC represents the onset of SAL episodes, which are 
associated with large warm anomalies located just off the west coast of Africa. The 
second PC has a dipole pattern that represents opposite phases of the evolution of the 
SAL. In the positive phase, the SAL tends to migrate southwestward into the subtropical 
to tropical Atlantic Ocean. In the negative phase, the SAL tends to migrate into the 
subtropical to mid-latitude Atlantic Ocean, the northwest Africa and southwest Europe. 
The phase of the second PC is related to the direction of transient eddy winds at the west 
coast of Africa around 15-20°N. In its positive phase, the transient eddy winds 
strengthen the climatological northeasterly, which blows offshore. In the negative phase, 
the eddy circulation pattern is reversed and the offshore flow at the west coast of Africa 
around 20°N is weakened. 

The NCEP meteorology is used to drive the CARMA model to simulate the 
corresponding dust distribution and deposition. We find that dust deposition is more 
likely to occur over northwest Africa and southwest Europe during the negative phase of 
the second PC than during the positive phase. The dust AOT anomalies simulated by 
CARMA have patterns in agreement with those seen in MODIS AOT anomalies, 
although the CARMA model underestimates the total AOT. Compared with the AIRS 
temperature retrievals in September 2002, NCEP temperatures at 850 hPa reproduce the 
migration of the warm anomalies associated with the SAL of our Episode IV. These 
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agreements between independent satellite measurements and model simulations increase 
our confidence in the reliability of our analysis and the quality of NCEP reanalysis and 
CARMA simulation. 

Several issues discussed in this paper need further research. We have identified that the 
evolution of SAL described by the second PC is associated with the strengthening of the 
easterlies at 700 hPa around 15-20°N, because of the thermal wind relation. Several 
studies [Carlson and Prospero, 1972; Cook, 1999; Dunion and Veldon, 2004; etc.] have 
indicated or implied the relation between the SAL and the African Easterly Jet, a 
component of the West African Monsoon (WAM). It will be interesting to further 
investigate the relationship between the evolution of the SAL and the subseasonal 
variation of the WAM. The traditional WAM index (WAMI) is defined by the difference 
of 200-hPa monthly mean zonal winds and the monthly mean modulus of the 850-hPa 
winds, averaged over 20°W-20°E and 3°-13°N [Fontaine et al., 1995; Moron et al., 2004; 
Palmer et al., 1992], Because we are studying the subseasonal scale variation of the 
circulation, we apply the same definition of the WAMI but replace monthly mean data 
with our daily noon-time wind anomalies. We have correlated our second PC time series 
with the daily WAMI and found significant correlations (-0.5 over 99% confidence level), 
if the PC time series lag behind the WAMI by 4-5 days. This implies a dynamical 
linkage between the variation of the WAM and the evolution of SAL. However, it is still 
not clear to us what the detailed physical mechanism behind these significant time-lag 
correlations is. It is necessary to further investigate whether the radiative effect of dust 
can influence the subseasonal variation of the monsoon-type wind patterns, and how the 
variation of the monsoon-type wind patterns affects the transport of dust. 

The second issue we would like to raise is about the dust radiative effect. It is clear from 
studies of simple models that heating from dust can alter the circulation pattern of the 
atmosphere [Stephens et al., 2004], Alpert et al. [1998] and Jones et al. [2004], by 
comparing the spatial distributions of temperature increments in data assimilations and of 
dust, suggested that dust has a warming effect at least over the eastern tropical North 
Atlantic Ocean. The correlation between the temperature increments and the dust 
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concentrations disappears in regions farther away from the African continent, possibly 
because of the lack of measurements for the assimilation increments. In our study, the 
CARMA dust is merely advected by the NCEP reanalysis winds and does not radiatively 
influence the winds. The NCEP winds are affected by dust radiative effects in the eastern 
Atlantic near the African coast through assimilation increment. However, there are less 
observational supports for the NCEP winds in the central Atlantic. The consistency of 
our results with satellite measurements even in the central Atlantic indicates that the 
relation between dust AOT anomalies and the temperature anomalies in regions farther 
away from the African continent is, perhaps, mainly determined by the dynamics. The 
warm anomalies originated from the African continent probably do not need the dust 
radiative heating to maintain it as it travels across the North Atlantic Ocean, and the 
warm anomalies and the dust anomalies may not be collocated with each other. The 
question of how far the radiative effect of dust seen in the eastern tropical Atlantic can 
extend and still be effective needs further research with a fully-coupled model. 
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Figure Captions 

Figure 1. Time series of the first (top) and the second (bottom) principal components of 
the variation of 850 hPa NCEP reanalysis temperature anomalies for June-September 
2002 over 10°-20°N and 20°-80°W. The abscissa is the Julian day number starting 
from June 1, 2002. The ordinate is scaled by the standard deviation of the time series 
(e.g., ±1 means at the ±l-a levels). The contributions of the PCs to the total 
variability are also shown in the figure titles. The arrows in the top panel show the 
four main episodes being discussed in the text. 

Figure 2. Geographical distributions of the first two (left panels) empirical orthogonal 
functions (units in K) of the variation of 850 hPa NCEP reanalysis temperature 
anomalies for June-September 2002 over 10°-20°N and 20°-80°W. The right panels 
show the CARMA dust AOT (at 0.55 pm) and the NCEP transient eddy winds at 700 
hPa projected onto the time series of the first (top right) and second (bottom right) 
PCs. 

Figure 3. Time-lag evolutions of the 4-day averages of the first (solid lines) and second 
(dashed lines) PC time series, starting from the onset dates of four main warm 
episodes. The ordinates are scaled by the standard deviations of the time series. 

Figure 4. The second-day lag composites of the 850-hPa NCEP temperature anomalies 
(in K) and 700-hPa transient eddy winds (left) and of CARMA dust AOT (right). The 
top panels show the composites for Episode I and the bottom panels show the 
composites for Episode II. 

Figure 5. Time-lag composites of the 850-hPa NCEP temperature anomalies (in K) and 
700-hPa transient eddy winds for Episode III. The zero-, second-, fourth-, and sixth- 
day lags are shown. 



22 


Figure 6. Time-lag composites of the CARMA dust AOT computed at 0.55pm for 
Episode III. The zero-, second-, fourth-, and sixth-day lags are shown. 

Figure 7. Similar to Figure 6 but for MODIS AOT at 0.55 pm. Missing data are set as 
negative and in white color. 

Figure 8. Similar to Figure 5 but for Episode IV. 

Figure 9. Comparisons of time-lag composites computed from the 850-hPa NCEP 
temperature anomalies and the 850-hPa AIRS temperature anomalies (units in K) for 
Episode IV. In this figure, the anomalies are computed as the difference of the 
temperature fields and the September mean. The top panels show the results for the 
second-day lag and the bottom panels show the results for the fourth-day lag. 

Figure 10. Similar to Figure 6 but for Episode IV. 

Figure 1 1 . Similar to Figure 7 but for Episode IV. 
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Table 1. Dust deposition over the eastern tropical North Atlantic Ocean (20°-40°W and 
0°-30°N) and the northwest Africa and southwest Europe (0°-20°W and 30°-45°N) 
for the four SAL episodes in June-September 2002. The first row in each episode is 
the average of dust deposition for the first four days (day 0 to 3) of the episode; the 
second row is the average over the next four days (day 4 to 7). Relative distributions 
over each region are also shown in percentage. 


Episode 

Z 1 I 

E. Atlantic (%) 

N.W. Africa and 
S. W. Europe 
(%) 

Averaged 
2 nd -PC Index 

I 

0.18 

34 

66 

-0.58 


0.19 

27 

73 

-1.69 

II 

0.19 

61 

39 

0.04 


0.31 

53 

47 

0.12 

III 

0.24 

85 

15 

1.55 


0.21 

50 

50 

-2.07 

IV 

0.27 

94 

6 

0.84 


0.19 

95 

5 

0.07 


















T850 PC -2 10.5% 


AOT and 700 hPa Winds on PC— 2 





Scoled Unit 













CARMA APT Episode 1 Lag: 2 


NCEP T850 Episode. 1... log: 2 


CARMA AOT Episode 2 Log: 2 


NCEP T850 















WsStea? 




NCEP T850 Episode 3 ..Lag: 0 


NCEP T850 Episode 3 Log: 


NCEP T850 Episode _ 3 Log: 4 


NCEP T850 Episode3 Log: 6 





MtiMii 


Swi 













NCEP T850 Episode, ft Log: 2 


NCEP T850 Episode, 4 Log: 0 


NCEP T850 Episode 4 











Figure 10 




MODS APT Episode 4 Lag: 0 


MODIS APT Episode 4 Lag: 2 


MODIS APT Episode 4 Lag: 4 


MODIS AOT Episode 4 Lag: 6 






